Universality of the Plastic Instability in Strained Amorphous Solids 
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By comparing the response to external strains in metallic glasses and in Lenard- Jones glasses we 
find a quantitative universality of the fundamental plastic instabilities in the athermal, quasistatic 
limit. Microscopically these two types of glasses are as different as one can imagine, the latter 
being determined by binary interactions, whereas the former by multiple interactions due to the 
effect of the electron gas that cannot be disregarded. In spite of this enormous difference the plastic 
instability is the same saddle-node bifurcation. As a result the statistics of stress and energy drops 
in the elasto-plastic steady state are universal, sharing the same system-size exponents. 



The nature of the fundamental plastic instability in 
simple models of amorphous solids with binary inter- 
actions was elucidated in recent years [l-'sj. The fun- 
damental plastic instability is best studied in athermal, 
quasi-static deformation of an amorphous solid such as 
to eliminate the effect of thermal fluctuations and finite 
strain rates. Contrary to some pictures suggesting that 
plastic instabilities occur at pre-determined 'weak' sites 
called 'shear transformation zones' H, [iqI? it was discov- 
ered that the plastic instability occurs as the vanishing of 
an eigenvalue of the Hessian matrix H where the latter 
is defined as 
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where U is the total potential energy which is a function 
of the positions of all the N particles in the system. Be- 
sides the three obvious zero eigenvalues associated with 
Goldstone modes, all the other eigenvalues of H are real 
positive as one expects from a real positive-definite sym- 
metric matrix. The plastic instabilities exhibit a univer- 
sal nature for all the studied binary glasses: at zero ex- 
ternal strain (7 = 0) all the eigenvalues of H besides the 
Goldstones are positive, and the low lying eigenvalues are 
associated with extended eigenfunctions. As the exter- 
nal strain is increased one eigenvalue begins to go down 
towards zero, and at the same time its associated eigen- 
function becomes localized. The eigenvalue, denoted Ap, 
hits zero at via a saddle node bifurcation (such that a 
minimum in the potential energy landscape hits a saddle) 
and the way that it does so is universal, i.e. 
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This scaling law has many interesting consequences. One 
immediate consequence is that the barrier AE between 
the minimum in which the system resides and the saddle 
that it eventually collides with tends to zero like [11 1 
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Less obvious (and maybe more interesting) is the conse- 
quence of these scaling laws on the system size depen- 
dence of the statistics of stress and energy drops in the 



elasto-plastic steady state. Once the stress level reaches 
the 'yielding transition' and the system settles in a steady 
state, every plastic instability is associated with a sta- 
tionary statistics for the drops Aa in the stress and AU 
in the potential energy. The long strain averages of these 
quantities scale with the system size according to [l, 0: • 



(AU) - A/-", (Act) - 



(4) 



It was shown that the numerical values of the exponents 
are determined by two ingredients. The first is the scal- 
ing laws ([2]) and (jSj. The second one is the assertion 
that after the yielding transition the probability to see a 
zero value of AE is not zero, in distinction from the un- 
strained solid in which this probability is zero, Il3| . 
The two ingredients combined result in a = 1/3 and 
/3 = —2/3 universally for all binary glasses and in both 2 
and 3 dimensions [l2|. In addition, it was shown that the 
statistics of plastic events continue to be dominated by 
the fundamental plastic instabilities discussed here also 
for finite strain rates and temperatures up to about 2/3 
of T^, see Ref. [14]. 

In this Letter we study the fundamental plastic insta- 
bility in a totally different class of glasses, i.e. metallic 
glasses where the microscopic interaction is very differ- 
ent from binary since the conducting electron cloud me- 
diates interactions beyond the pair-wise. We study the 
response of a model of the metallic glass Cu^qZt^^ in 2 
and 3 dimensions using the Embedded Atom Model, and 
show quantitative universality with the binary glasses. 
In spite of the tremendous difference in microscopic in- 
teractions, the scaling laws (|4]) will be shown to repeat 
verbatim, with the very same scaling exponents. 

In order to simulate the consequences of the electron 
cloud we use a many-body potential. This potential rep- 
resents the cohesive energy of an atom by the local elec- 
tron density at that atom, the latter quantity itself being 
determined by the neighbors of the atom. A pure pair- 
wise interaction is added in order to represent the elec- 
trostatic repulsion (see pjj and references therein for a 
more detailed description). The expression for the total 
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FIG. 1: Upper panel: typical stress vs strain curve for a sys- 
tem of 4000 particles. Lower panel: a typical potential energy 
vs strain curve for the same system as in the upper panel. 



potential energy is 



' mn — ' n 'ml 
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Here ^ is the pair- wise potential, F (pm) is the many- 
body term representing the energy of the mth atom and 
Pm is a measure of the local electron density at the lo- 
cation of the mth atom, due to contributions from the 
neighboring atoms [HI. The functional form of F {pm) 
is taken as customary in the literature, cf. lH, [l7| i.e. 
the semi-empirical approximation F {pm) ~ —^/p^^ We 
employed the functional forms of (j) and ip from Ref. [HI 
where fits were obtained for this model to density func- 
tional quantum mechanical calculations, for the case of 
the bulk metallic glass Cu46Zr54. For the sake of compu- 
tational speed we cut-off the functions ^ and ^ to go to 
zero smoothly (up to second derivative) at a finite dis- 
tance. Explicitly, we used the following functional forms 
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FIG. 2: Color online: Potential Energy per particle versus 
strain for systems of different sizes between N= 400, and N — 
8000. Note the energy peak and then dip before the elasto- 
plastic steady state is established, this seems particular to the 
metallic glass and not commonly seen in binary glasses. 
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Here the coefficients A, B and C have been added to 
achieve smooth first and second derivatives of (j) and tjj 
at the cut-off. The values of the parameters q a and 
c in equations ([7]) and ([8]) are taken from [18]. Note that 
cutting off the binary functions at a finite distance does 
not remove the multi-body interaction because the force 
between any pair or particles is configuration dependent 
through the second term in the total energy. For our 
purposes it suffices to smooth the potential up to second 
derivative since we are interested in the Hessian matrix. 
We use a (scaled) value of cut-off 1.707 corresponding 
to an absolute value of 4.5A used in [l8|. The scaled 
value of density is 0.951 obtained from the data in [18] 
corresponding to 2000 atoms occupying a volume of 38.5 
nm^ (close to Tg in [HI). All the simulations reported 
here are in three dimensions, keeping a constant density 
in a 3-dimensional box with periodic boundary conditions 
in all directions. 

We have put our model metallic glass under a shear 
strain using the usual algorithm for an athermal quasi- 
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FIG. 3: Color online: The scaling of (At/) (upper panel) and 
(Act) (lower panel) with the system size N . Symbols indicate 
numerical data while the dashed line is a fit. The dashed 
curves are power laws with exponents a — 0.330 and /3 = 
—0.647 respectively, in good agreement with the power laws 
dll) which predict universal exponents a = 1/3 and /3 = —2/3. 

static process [7]. Fig. [T] demonstrates a typical stress 
vs. strain and potential energy vs. strain curves for a 
metallic glass consisting of 4000 particles. We see the 
usual linear regime in the upper panel where the stress is 
linear in the strain, and the corresponding regime where 
the potential energy is quadratic in the strain in the lower 
panel. Every drop in the curves represents a plastic fail- 
ure, and we are interested in the statistics and system 
size dependence of those. Indeed, in Fig. [2] we show the 
potential energy per 'particle vs strain for systems of dif- 
ferent sizes from = 400 to TV = 8000. We see the 
decline in the amplitude of the fluctuations but also a 
small reduction in the average energy in the elasto-plastic 
steady state. Note the energy peak and then dip before 
the elasto-plastic steady state is established, this seems 
particular to the metallic glass and not commonly seen 
in binary glasses. 

The scaling laws with respect to the system size are 
examined in Figs. [3l We see that Eqs. |4]are reproduced 
by the data with the scaling exponents being very close 
to the theoretical values a = 1/3 and /3 = —2/3. This 
demonstration of universality is determined by the fun- 
damental plastic instability which, in this case of metallic 
glass, in spite of the huge difference in microscopic inter- 
actions, is the same saddle node bifurcation as in all the 



FIG. 4: Color online: The lowest eigenvalue of the Hessian 
matrix as it drops to zero at a typical plastic event in the 
elasto-plastic steady state. The agreement with the saddle- 
node bifurcation scenario (dashed red curve) seems perfect. 

simple cases of binary glasses. To show that this is so 
we need to diagonalize the Hessian whose computation 
is somewhat delicate in the present case. The analytic 
form of the Hessian for our Embedded Atom Model is 
a bit cumbersome to write here explicitly, but it can be 
found in A typical calculation of the lowest eigen- 
value for the case N = 400 is shown in Fig. HI We see 
how the eigenvalue goes to zero at a value 7p precisely 
as expected for a saddle node bifurcation, i.e. according 
to Eq. (j2j). 

The conclusion of this study is two- fold. First, we learn 
that the nature of the precise microscopic model does not 
matter much in determining the fundamental plastic in- 
stability of amorphous solids. Once coarse-grained, or 
up-scaled to provide us with an energy landscape, the 
plastic instability occurs when the strain distorts the 
landscape such that a local minimum collides with a near- 
by saddle. Such a collision is a generic saddle node bi- 
furcation that will be seen as an eigenvalue of the Hes- 
sian matrix nearing zero with a square-root singularity. 
The power of genericity of bifurcations exceeds the de- 
tails of the inter-particle potential. Secondly, we learn 
once more that it is futile to seek 'shear transformation 
zones' in the material. At rest, unstrained, the mate- 
rial does not exhibit any pre-existing local weak points 
where the plastic event will take place. The plastic event 
is a result of the straining process, in which an extended 
eigenfunction of the Hessian matrix gets localized at the 
position of the plastic event. This lesson is as true for 
simple glasses as for the metallic glass. We conclude this 
letter by showing the eigenfunction associated with the 
lowest eigenvalue at 7 = and at 7p in a 2-dimensional 
metallic glass (for ease of display, cf. Fig. [5]). In equilib- 
rium the eigenfunction is perfectly extended, and there is 
nothing in it that can predict the position of the incipient 
instability. Once localized, the eigenfunction exhibits the 
typical quadrupolar structure tha t is associated with the 
well-known Eshelby inclusion [20], which is universal for 
elastic materials independent of the microscopic details. 
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FIG. 5: Color online: A 2- dimensional demonstration of the 
nature of the plastic instability. Upper panel: at 7 = the 
eigenfunction associated with the smallest non-zero eigenvalue 
is perfectly extended without any sign of the position of the 
incipient instability. Lower Panel: the localized eigenfunction 
associated with the very same eigenvalue when it approaches 
zero, (see Fig. SI). 
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